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Abstract. The statistical properties of heart beat intervals of 130 long-term surface 
electrocardiogram recordings during atrial fibrillation (AF) are investigated. We find 
that the distribution of interbeat intervals exhibits a characteristic exponential tail, 
which is absent during sinus rhythm, as tested in a corresponding control study 
with 72 healthy persons. The rate 7 of the exponential decay lies in the range 3- 
12 Hz and shows diurnal variations. It equals, up to statistical uncertainties, the 
level of the previously uncovered white noise part in the power spectrum, which 
is also characteristic for AF. The overall statistical features can be described by 
decomposing the intervals into two statistically independent times, where the first 
one is associated with a correlated process with 1// noise characteristics, while the 
second one belongs to an uncorrelated process and is responsible for the exponential 
tail. It is suggested to use 7 as a further parameter for a better classification of AF 
and for the medical diagnosis. The relevance of the findings with respect to a general 
understanding of AF is pointed out. 
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1. Introduction 

Atrial fibrillation (AF) is the most common arrhythmia of the heart and 
leads to an impairment of physical ability and to a reduction of quality 
of life (for a review, see e.g. [1]). Connected with these degradations is 
an increased risk of thromboembolic complications, of strokes [2] and 
of mortality [3]. The prevalence of AF in the industrialised countries 
currently is about 1% with respect to the total population and gives rise 
to tremendous costs for diagnosis and therapy. The risk to develop AF 
increases with age: For persons older than 40 years 2.3% are affected, 
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for persons older than 60 years 5.9%, and for persons older than 80 
years 8% [4]. Because of the rise of life expectancy in the industrialised 
countries, the number of persons with AF is likely to increase further. 

To support the choice between currently employed therapeutical 
measures, such as pharmaceutical treatment with antiarrhythmica, AV- 
nodal ablation or modulation, pulmonary vein or isthmus ablation, and 
cardioversion, it is necessary to develop a good classification of AF. 
The size of the left atrium is to date the main parameter to support 
diagnosis and decisions for treatment options [5, 6], although its role 
in predicting outcome after a cardioversion is discussed controversially 
[7, 8]. Further classifications are based on the averaging of P waves 
[9, 10], on the amplitude of the fibrillatory waves in the surface electro- 
cardiogram (ECG) [11], or on the morphology of the intra-atrial EGG 
[12]. The number and degree of spatial uniformity of wave fronts was 
used for a classification into 3 different AF types [13]. In recent studies 
short cycle lengths of fibrillatory waves were found to correlate with 
the risk, that paroxysmal AF changes into persistent AF [14, 15] and 
that arrhythmica fail to induce conversion to the sinus rhythm [16, 17]. 

All these established methods use information from intra-atrial sig- 
nals (or atrial contributions in the surface ECG) to classify AF. How- 
ever, also ventricular beats show characteristic signatures of AF [18, 
19, 20], parts of them have been suggested for the detection of AF 
[18, 22]. In this work we will show that the distribution of RR intervals 
during AF exhibits a characteristic exponential tail with a decay rate 
7 in the range 3-12 Hz. We show that this rate can be determined also 
in the power spectrum. Our analysis suggests to use 7 as an additional 
classification parameter for AF. 



2. Characteristic features of atrial fibrillation 
in RR interval distributions 

Our study is based on RR sequences of 130 AF patients, which 

typically comprise 10^ beats (corresponding to 24h EGG recordings 
[?]). The probability density p(t) of RR intervals for one representative 
patient is shown in Figs. la,b. We find that for large r, the decay of 
p{t) can be well fitted to an exponential, 

p{r) PooCw{-1t) , (1) 

with decay rate 7 and amplitude factor p^o- This behaviour becomes 
particularly clear in the semi-logarithmic representation in Fig. lb. 

We find that eq. (1) is a generic feature for AF patients. To quantify 
the significance of (1) we perform a Kolmogorov-Smirnov test [23] for 
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the RR intervals in the tail of p{t). For sufficient statistical meaning, 
wc require the tail region to encompass at least 0.3 s and at least 2% of 
the RR intervals. We find that 106 of the 130 AF patients (81.5%) pass 
the Kolmogorov-Smirnov test with a standard significance level of 10%. 
Even better values are obtained by a more sophisticated procedure for 
identifying the parameters 7 and Poo from an optimal fitting region for 
eq. (1) as described below. 

One can expect that eq. (1) is also specific for atrial fibrillation, since 
(1) it is not known to be a particular feature of healthy persons. To 
perform a countercheck, we analyse RR sequences of 72 healthy persons, 
which were taken from the PhysioNet database [24] . The result for one 
representative person is displayed in Figs. lc,d. Clearly, no exponential 
tail can be identified. Only 6 of the 72 healthy persons (8.3%) would be 
accepted according to the Kolmogorov-Smirnov test as specified above. 

The rough method to identify the tail yields approximate 7 values 
in the range 3-12 Hz. However, to determine the parameters poo and 
7 for each patient more precisely, a "best fitting" interval [ti , T2\ is 
calculated for each patient, in which the Kolmogorov-Smirnov deviation 
is smallest. In addition to the requirement that the tail region should 
encompass at least 2% of the total number of RR intervals, we vary both 
Ti and T2 under the constraints (i) T2 — T1 > 0.3s and (ii) ti > 0.7s. The 
first constraint (i) guarantees that, in the case of small 7, the width 
of the fitting region covers about one decay time I/7. In the case of 
large 7, it ensures that the fitting interval covers several I/7 so that it 
is not influenced by short-time fluctuations. The second constraint (ii) 
in combination with (i) ensures that the fitting region lies right to the 
maximum of p{t). In the best fitting region, 125 of the 130 AF patients 
(96.2%) pass the Kolmogorov-Smirnov test with significance level 10% 
and these are taken into account in the following analysis. The fitting 
region has a mean width of 0.42 s and on average comprises 8.6% of the 
RR intervals. The histogram of 7 values for these patients is displayed 
in Fig. 2. The values lie in the range 2-12 Hz and have a mean 5.4 Hz. 

While almost all p{t) of the AF patients show the exponential tail, 
the shapes of the distributions differ outside the tail region. The shapes 
can be grouped into classes according to the number of local maxima 
in p{t) (unimodal or multimodal distributions). The physiological ori- 
gin of multimodal distributions is still under investigation and mostly 
associated with the presence of different conduction pathways through 
the AV node [25, 26, 27]. 

A full account for the unimodal distributions can be obtained by 
assuming that r results from a superposition of two statistically in- 
dependent times, r = t' + rj, where r' is drawn from a Gaussian 
distribution i/j{t') with mean tq and variance ctq, and r] is drawn 
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from an exponential distribution 0(ry) = 7exp(— 77?) [28]. In order for 
the times r' to be positive one could formally introduce a truncated 
Gaussian, but because it turns out that iP{t') is sufficiently sharply 
peaked at tq, we can ignore the truncation. Then we obtain 



p(r) = r dr' i;{T') (j){T - t') 

J — OD 



= 7 explTTG + 2) ^^^^ ^ tq + 7^ - r ^ 

with mean r = tg + 7^^ and variance cr^ = ctq + 7^^. Hence, given the 
7 values as determined above, fits to eq. (2) can be easily performed 
by just calculating f and a^. These fits are in good agreement for all 
48 AF patients with a unimodal distribution. For the representative 
example in Figs. la,b. this is demonstrated by the solid lines. For large 
RR intervals the distribution follows eq. (1) with 

f / 7^2 3\ 
Poo = 7 exp 7TG H = 7 exp 7r + — (3) 



As shown in Fig. 3 the Poo determined by the linear regression follow 
closely the theoretical curve predicted by eq. (3). This allows us to scale 
the exponential tails for these AF patients onto a common master curve, 
see the inset of Fig. 3. The tg and ctq lie in the ranges 0.3 — 0.85 s and 
0.04 — 0.17 s respectively. 

The decay rate 7 can be viewed as a new parameter for a further 
classification of AF. As such, it may be subjected to diurnal variations 
as, for example, the fibrillation rate [29]. To test this, we next study 
correlations between 7 and tg, which can be expected to be larger 
during night than day time. The sequences of RR intervals are split 
into segments of 2^'^ intervals (corresponding to typically 4.5 hours). 
For good temporal resolution, the analysis is performed on overlapping 
segments shifted by 2^^ beats (corresponding to typically 1 hour). In 
each segment i, ji is determined from the RR interval distribution 
belonging to this segment as described above, and rG,i = fi — is 
calculated from the the mean fj in the segment. 

Figure 4a shows the results for a representative patient. Clearly, 
time periods with larger tq correspond to periods of smaller 7 and vice 
versa. To quantify this behaviour we calculated the cross correlation 
coefficient C for all AF patients. The histogram of C values in Fig. 4b 
reflects pronounced anti-correlations. The mean correlation coefficient 
isC = -0.86. 
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3. Relation to properties of the power spectrum 

Properties of the power spectrum S{f) of the RR interbeat intervals 
[30] during AF have been studied previously by some of the authors 
and coworkers [20, 21]. For low frequencies S{f) exhibits a 1// type 
behaviour as known also for healthy persons. At higher frequencies, 
however, striking differences between the power spectra of healthy per- 
sons and AF patients arc found (see Fig. 5a for a representative ex- 
ample). While the 1// type behaviour continues up to the highest 
frequencies for the healthy persons, there is a crossover to an extended 
part with a behaviour close to white noise at high frequencies for the 
AF patients. The question arises if the two characteristic features, i.e. 
the exponential tail in the distribution of RR intervals and the white 
noise part in the power spectrum, are interrelated due to a generic 
underlying mechanism. 

Motivated by the successful decomposition of the RR intervals r into 
two contributions r' and r/, we suppose the exponentially distributed 
times 1] to be uncorrclatcd and hence associated with the white noise 
part in S{f). By contrast the time r' should reflect a strongly correlated 
process associated with the 1/f part. Accordingly we write S(f) = 
Sr'if) + Srfif), where S^'if) ~ Sof~°' and S'^(/) ~ are the power 
spectra of the two contributions. Then the crossover frequency is given 
by /x = (So/Soo)^/". The strength S^o of the white noise should be 
given by the decay rate 7 (since = 'S',,(/) = Sn{f 0), and the 
zero-frequency limit of a power spectrum equals the variance), 

5oo = l/7'- (4) 

— 1/2 

Using (4), we determined 7pow = by averaging Sri{f) over the 

modes belonging to the 100 largest frequencies [30]. The result is shown 
in Fig. 5b in comparison with the decay rate 7 obtained from p{t). A 
fair agreement is found for all 48 AF patients. Also, 7pow follows closely 
the diurnal variations of 7. This can be seen by the dashed line and 
the hatched bars in Figs. 4b and 4c, respectively. Nevertheless, when 
doing the time-resolved analysis, the 7 obtained from the RR interval 
distribution and the power spectrum are not as close as in Fig. 4b. 
This is due to the fact that a precise determination of 7 from the RR 
interval distribution becomes rather difficult for a comparatively small 
number of RR intervals. By contrast, the white noise part in the power 
spectrum can be clearly identified already for segments comprising 2^° 
beats only. As a consequence, the power spectrum should be preferred 
for the determination of the new parameter 7 in medical applications. 
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4. Conclusions 

In summary we have studied the distribution and correlation properties 
of heart beat intervals based on long-term (24h) EGG recordings during 
AF. A characteristic exponential tail in the distribution of heart beat in- 
tervals was found, which is absent in healthy controls and can therefore 
be used as a novel way to characterise AF. Wc quantify the tail with 
the decay rate 7 and find values in the range 3-12 Hz. Surprisingly, 
these 7 values can also be determined from the white noise part of 
the power spectrum. The findings can be described by decomposing 
the RR intervals into two statistically independent times, where the 
first one belongs to a correlated process with the typical 1// noise 
characteristics (as found also for healthy subjects), while the second 
one belongs to an uncorrelated process. The uncorrelated process is 
responsible for the occurrence of the exponential distribution and the 
white noise. Over the day the 7 values vary slowly and these variations 
are anticorrelated with the mean interval belonging to the correlated 
process. While the time-resolved analysis is rather difficult to perform 
based on the distribution function, it is straightforwardly carried out 
based on the power spectrum. 

It is remarkable that the range of 7 values as well as its diurnal 
variations are comparable to the range and diurnal variations of the 
fibrillation rates typically found in medical studies [31, 32, 29]. While 
it is currently unclear if this is a mere accident or if there exists a 
deeper connection, the potential of the new parameter 7 for the clinical 
diagnosis should be evaluated. This parameter is readily available from 
the surface EGG and accordingly easier to obtain than data from intra- 
atrial recordings. Similar to the fibrillation rate, 7 may provide a useful 
measure to estimate the outcome of a cardioversion or of pharmaceu- 
tical treatments. Among others, it could also help to predict the risk 
that patients with paroxysmal AF develop chronical AF. We hope that 
our findings will stimulate further research in this direction. 

Apart from this issue of practical importance, the results require an 
explanation based on physiological mechanisms connected with AF. For 
such explanation both the statistics of atrial impulses and the trans- 
duction of them through the atrio- ventricular node are expected to play 
a key role. On may assume that the observed statistics directly reflect 
statistical features of the intra-atrial beats. However, typical recordings 
in the right atria show a rather narrow distribution of the intervals 
between successive atrial impulses and therefore do not support this 
idea. On the other hand, it has been suggested that blocked atrial 
impulses can prolongate the (absolute or relative) refractory period of 
the atrio-ventricular node, and this process could be a possible reason 
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for the peculiar statistics of the long heart beat intervals during AF. 
Further studies, however, are needed to clarify this possibility. 
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Figure 1. Probability density p(r) of RR intervals r for one representative AF 
patient on (a) linear and (b) semi- logarithmic scale. In the calculation a fixed bin 
width of 0.03 s was used. The lines mark the fit to eq. (2) with the exponential decay 
for large r (cf. eq. (1)). Corresponding plots on (c) linear and (d) semi-logarithmic 
scale of one representative healthy person show no signature of an exponential decay. 
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Figure 2. Histogram of 7 values with mean 5.4 Hz and standard deviation 2.1 Hz 
for all 125 patients. 




Figure 3. Semi-logarithmic plot of poo/7 for the subgroup of 48 patients as a function 
of 7f +7^o"^/2. The straight line marks the theoretical behaviour according to eq. (3). 
The inset shows the rescaled probability density for this group of patients. 
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Figure 4- (o-) Time dependence of tq and decay rate 7 in segments of 2^ beats 
during AF (solid line: 7 calculated from the RR interval distribution; dashed line: 
7 = 7pow calculated from the power spectrum, see Sec. 3). The segments were shifted 
by 2^^ beats, and rc.i and 7^ are assigned to the starting point of each segment i. (h) 
Histogram of cross-correlation coefficients C — ~ ''^)(7» ~ 7)/(-^^7^t), 

where A'' is the number of segments, f — X)i!=i '''G,i/N, 7 = Ji/N, and A,-, 

A-y are the corresponding standard deviations (full bars: 7^ calculated from the 
RR interval distribution; hatched bars: 7^ calculated from the power spectrum, see 
Sec. 3). 
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Figure 5. (a) Power spectrum of RR intervals r for one representative AF patient on 
a double-logarithmic scale. The spectrum exhibits a distinct crossover at a frequency 
/crit — (100 beats) For frequencies / < /ciit the spectrums exhibit the typical 
1/f behaviour and a white noise characteristic for / > /crit- (b) "/ vs. 7pow for the 
group of 48 patients with unimodal p(t). The red line marks the theoretical relation 
7 = 7pow. 
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